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I. INTRODUCTION 


Quantum graphs have proven to be useful models in various areas of mathematical and 
theoretical physics. From a mathematical point of view, they combine the simplicity of a 
(quasi) one-dimensional system with the complexity of a graph-like structure. It is exactly 
this underlying complexity that turned quantum graphs into particular important models 
in the held of quantum chaos by showing that eigenvalue correlations in quantum graphs 
are generically the same as in quantum systems with chaotic classical limit^'^. Intuitively, 
the chaotic behavior stems from the scattering of the quantum particle in the vertices of the 
graph. Therefore, in order to provide a useful model for the investigation of many-particle 
quantum chaos, many-particle systems on (hnite, compact) quantum graphs with localised 


many-particle interactions were introduced in Refs. Q and 0. Localised in this context means 
that the particles interact only in the vicinity of the vertices of the graph, i.e., one particle has 
to sit at a vertex whereas the other particles are close to it. The important consequence of 
such singular many-particle interactions is that the scattering of the particles and hence the 
dynamics of the system are affected in a way encoding genuine many-particle correlations. 

Also, apart from the fact that the model which will be discussed in this paper originated 
from the held of quantum chaos, it is also an interesting model in its own right. From the 
point of view of ^plications, singular many-particle interactions on graphs were already 
discussed in Ref. l5| in order to investigate the effects of short-range many-body interactions 
on, e.g., the conductivity of nanoelectronic devices. More explicitly, our model can be 
understood as a toy-model to investigate a system of two electrons moving in a wire which 
is normal-conducting except for a relatively small part where it is super-conducting^>^. In 
the super-conducting part the electrons are then interacting due to the pairing effect in 
superconductivity (Cooper pairs). 

From a more theoretical point of view, it is well-known that quantum many-body prob¬ 
lems are generally hard to solve and, indeed, there are only few models which are explicitly 
solvable^. One important model in this respect is the Lieb-Liniger model (originally formu¬ 
lated for an A^-particle system), see Ref. 0, whose Hamiltonian, in the case of two particles. 


is formally given by 


H = 




+ a6{x — y) , 


( 1 ) 


dx‘2 dy"^ 

where a G M+ is the interaction strength and 6{x) the Dirac delta function. Explicitly 
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solvable for this model means that its eigenfunctions can be (if it is considered on an interval 
together with periodic boundary conditions) explicitly given and the spectrum can, at least 
in principle, be calculated. We note that versions of the Lieb-Liniger model on graphs were 


discussed in Refs, y and 


10 


in the case of star graphs and in Ref. [2 in the case of general 


compact graphs. We will see below that the formal Hamiltonian of our model resembles the 
Hamiltonian ([1]) to some extent. The difference is that the two-particle interactions which we 
will consider are even more singular in the sense that they do not solely depend on the relative 
position of the two particles. It is important to realise that, although the Lieb-Liniger model 
was originally considered as being of theoretical interest only, it is nowadays recognised to 
describe realistic gases in one dimension very well at low temperatures^^. Hence, in the same 
spirit as for the Lieb-Liniger model, we consider the model to be discussed as an interesting 
model for theoretical as well as for practical purposes. 

This paper is organized as follows; In Section |TT] we introduce our model and give a 
precise mathematical formulation, i.e., we establish a quadratic form which is then shown 
to characterise a self-adjoint operator (the Hamiltonian of our system) uniquely. Section HTTI 
is devoted to the spectral analysis of this operator, i.e., we will characterise its essential 
spectrum and establish conditions for which there exists at least one eigenstate below the 
essential spectrum. We will also provide upper bounds on the number of such eigenstates. 
Finally, in Section IIVI we will establish (pointwise) upper bounds on the ground state 
eigenfunction proving an exponential decay as -f —)■ oo. 

II. THE MODEL AND ITS HAMILTONIAN 


The model we would like to discuss in this paper consists of two (distinguishable) interact¬ 
ing particles where the one-particle conhguration space is given by the half-line M+ = (0, cxd). 
The two-particle interactions shall be of singular type, i.e., the formal Hamiltonian of the 
system reads 

02 

^ ^ ^ ^ ^ 

where v{x,y) = v{y,x) is some symmetric interaction potential. On a physical level we can 
read off from (E]) that the two particles interact only when at least one particle is situated at 
the origin. Furthermore, choosing v{x,y) such that suppn C B^{0) OM^, where B^{0) C 
is the open ball with radius £ > 0 around 0 G M^, we see that the particles interact only 


3 






whenever one particle is situated at the origin and the other particle is contained in a small 
neighborhood around it. 

On a mathematical level, the methods employed in Ref. Q show how the Hamiltonian 
(1^ can be realised via a quadratic form on Most interestingly, this quadratic 

form then corresponds to a variational formulation of a boundary-value problem for the 
two-dimensional Laplacian 

<92 


dx"^ dy"^ 

with coordinate dependent Robin boundary conditions. Indeed, defining (j{y) := —n(0,2/), 
the boundary conditions read 


|^(0,|/)+ a( 2 /)y?( 0 , 2 /) = 0 , and 


dip 

dn 


{y,0) + a{y)ip{y,0) = 0 , 


(4) 


where ^ denotes the inward normal derivative along 


Remark II. 1. We require < 7 ( 2 /) to be a measurable, essentially bounded function throughout 
this paper. 

The associated quadratic form is given by 


q[ip] = / dx 


(^{y) Wbv{y)? dy , 


(5) 


defined on Vq = i.e., the Sobolev space of order one. The (pbviy) are the boundary 

values of 99 G which are, according to the trace theorem for Sobolev functions (see 


Ref. 


12|), well defined. It is then readily verified that the form ([5]) is meaningful and we can 


establish the following result. 


Theorem II. 2. The form 


is closed and semi-hounded. 


Proof. For the proof we note that the restrictions of all functions in onto are 

dense in (see, e.g.. Ref. Il2h . This allows us to transfer trace estimates for bounded 


(Lipschitz) domains to the unbounded case, i.e., to R^. More explicitly, in the proof of 
Theorem 3.1 in Ref. Q, an estimate is used which extends to our case: we have 


1111x2(98^) — 4 ( ^ I|¥’IIl2(R^) + '^II^V^IIl2(r2 


( 6 ) 
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for all 5 > 0. We hence obtain 


(y{y) Wbv{y)? <^y 


< 4 ||cr(i/)||oo tII‘^IIl2(k2) + ^||V(/?||^2(k2 


(7) 


which implies 


q[^] > (1-45||a(i/)||oo) ||V(/?||^2 (r2; 


8|k(i/)llc 


L2(]R2 ) • 


( 8 ) 


Hence, choosing 6 small enough, we see that the form q[-] is bounded from below. 

Furthermore, the estimate ([7]) also implies that the form-norm is equivalent to the H^- 
norm. Since is complete, we conclude that q[-] is indeed closed. □ 

Remark II. 3. From equation (0) one directly concludes that 


inf (t(—A g-) > —32||(j 


2 

OO 


(9) 


by setting 6 = 


4|k(3/)l|c 


-. Note that a sharper bound will be derived later in Theorem IIII. 71 


Due to the representation theorem of quadratic forms (see, e.g.. Ref. Il3l) . Theorem III. 21 
implies the existence of a unique self-adjoint operator being associated to the form q[-]. This 
operator is the Hamiltonian of our system which we denote by —A^. 

Remark II. 4. We note that the case <7 = 0 corresponds to the so called Neumann-Laplacian 
on being self-adjoint on the domain := {yj G \ = 0 on 


III. SPECTRAL PROPERTIES OF -A^ 

In this section we characterise the spectrum of — Ao-, i.e., we describe the essential part as 
well as the discrete part of the spectrum. We will establish conditions on < 7 ( 2 /) for which there 
exists at least one eigenstate below the essential spectrum. Furthermore, we provide upper 
bounds on the number of eigenstates below the essential spectrum and specify conditions 
on (j{y) for which there is exactly one such eigenstate. 


A. On the essential spectrum 

In our first result we prove that the interval [0, cxo) is, in the most general case of (essen¬ 
tially) bounded < 7 ( 2 /), always contained in the essential spectrum. Furthermore, we prove that 
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the essential spectrum has no negative part if either u{y) becomes eventually non-positive 
or if |cr(?/)| converges to zero. 

Theorem III.l. For a{y) G one has 

[0, cx)) C aess(-A^) . (10) 

Furthermore, if either 

(i) there exists L > 0 such that u{y) < 0 for a.e. y > L , 


or 

(a) esslimj/^oo W{y) \ = 0 , (essential limit) 
then 

aess{-^a) = [0,Oo) . ( 11 ) 

Proof. To prove (fTOj) we employ Weyl’s characterization of the essential spectrum in the 
sense of quadratic forms^^. This means that, in order to prove that A > 0 is in the essential 
spectrum, we construct a sequence (</3ri)neN Tf^(M^) such that ||</?n|| = 1, 0 (weak 

convergence) and 

sup I |s((p„, U) - X{^n, m)l2(R^) 

as n —)■ oo and where s(-, •) is the sesquilinear form associated with ([S]). To construct such 
a sequence we consider the rectangle := [0, L„] x [j/„, yn + Bn] with some constants 
Ln,Bn,yn > 0 and choose (pn to be the normalized ground state eigenfunction of the two- 
dimensional Laplacian (with Dirichlet boundary conditions) on Rn, extended by zero to all 
of A direct calculation then gives 

/ TT^ TT^ A 

s{p>n,U) = -I- — 1 {ipn, U) . ( 13 ) 

2 

Hence, letting ^ oo, ^ A and ?/„—)■ oo fast enough, we see that ((pri)neN is an 

-^n 

appropriate sequence showing that A G (Tess(—Ao-). Hence [0, cx)) C cTess(—Ao-). 

To prove the second statement we use a bracketing argument. For this, we split 
into two disjoint subsets: we write = Dji U where := \ Dji and Dji = 

(0, R) X (0, R) with R> 0 some constant. As a comparison operator we consider the direct 
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sum © —Aq_ of Laplacians, the index N referring to (additional) Neumann boundary 

conditions along the inner two line segments of Du. Note that both operators, —A^^ as 
well as —A^_, can be dehned via their associated quadratic forms which are similar to ([5j) . 
More precisely, —A^^ is the unique self-adjoint operator being associated with 


?ib] := / \V(p\'^dx-[ a{y) \ipbv{y)\‘^ dy , (14) 

JDr JdM.%nBR{0) 

dehned on Furthermore, —A^_ is the unique self-adjoint operator being associated 

with 

gsM := / \Vip\‘^dx-[ a{y) \ipbv{y)\‘^ dy , (15) 

JU- JdRl\^R(0) 


dehned on 

A standard bracketing argument of operators, see eg. Ref. 


13 


, then implies 




(16) 


Since (Tess(—A^^) = 0 (note that — An is dehned on a bounded Lipschitz domain and hence 
has only discrete spectrum, see Ref. 


12l) we have 
infcTess( Aj^_) ^ mfcTess( A(j) . 


(17) 


Assume now that cr{y) fulhls condition [i): if we choose i? > L in the previous construction, 
—Aq_ is a positive operator since a{y) < 0 for a.e. y > R. Furthermore, using the same 
methods as in the proof of the hrst statement, we conclude that [0, cx)) C (Jess(—A(^_). 
Accordingly, by flTTj) and by flTOl) . one concludes that (Jess(—A^-) = [0, cx)). 

Finally, assume that (^{y) fulhls condition (ii): choosing again R > 0 in the previous 
construction large enough such that ||c’'(j/)||icx)(gR^\ 5 ^(o)) < £ one has (in the same way as 
outlined for — A^. in Remark III.3j) 

inf a(-A^_)>-32^2 . (18) 

Since e > 0 is arbitrary we conclude the statement while taking (ITT)) into account. □ 

Remark III.2. Regarding Theorem IIII.ll it is instructive to consider a case where a{y) is 
(essentially) bounded and where the essential spectrum has indeed a negative part. For this 
we choose a{y) = a > 0 which means that the associated quadratic form is given by 

qa[^]:= [ \V(p\‘^dx-a[ {(fbviy)]"^ dy , (19) 

./r2 
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being defined on Vq = The important observation is that the Laplacian — Aq- which 

is associated with go-[-] can be decomposed as 


A,, = -Ai.1) (g) 1 + 1 0 -Al^) 


( 20 ) 


where — A^-^^ is the (self-adjoint) one-dimensional Laplacian — defined on Do- = {(p G 
.ff^(M+) I (p'(0) -|- cr(p(0) = 0}. This follows from the fact that the quadratic form associated 

is 

■= f dx - a|(p(0)p , (21) 


d2 


with —aI^^ 


being defined on Also, it is well-known that — A^^ has exactly one eigenstate 

of negative energy at —cr^ with multiplicity one (see, e.g.. Ref. Q) and essential spectrum 
(Jess(—A ct^) = [0 ,cxo). Due to (l20|) one can employ standard results (see, e.g.. Refs. 17,18) 
to conclude that (Trf(—A^-) = {—2a‘^}, where —2cr^ is an eigenvalue of multiplicity one, and 
(Jess(—Ao-) = [—cr^,cxD). Furthermore, the (normalised) ground state associated with the 
eigenvalue —2(T^ is given by ip{x,y) = 


B. On the discrete spectrum 


In a second step we prove the existence of an eigenstate at negative energy whenever the 
two-particle interactions are attractive on average, assuming in addition that <j{y) is such 
that inf (Jess{—^a) = 0 (see condition (i) or {ii) as formulated in Theorem IIII.ll) . Attractive 
on average shall mean u{y) G L^(M+) and 

a{y) d?/ > 0 . (22) 

From a physical point of view, condition fl22p implies that the potential <j{y) can exhibit 
some oscillatory behavior while still maintains an eigenstate at negative energy. This, 
on the other hand, allows one to consider “more natural” interactions since repulsive inter¬ 
actions are always expected if the particles are too close to each other (hardcore repulsion). 



Theorem III. 3. Assume that aly) G L°°(M+) is such that inf cress(—Ag-) = 0 (compare 
with Theorem \III.l\) . If furthermore a{y) G L^(M+) and a(y)dy > 0 then there exists an 
eigenstate state below the essential spectrum. 











Proof. We use a test-function argument and consider the (unnormalised) function 


ipjr) = e , £ > 0 , 


(23) 


using polar coordinates. We calculate 


TT 


(l[Pe] = 2 Weir)? r dr - 2 a{y)e dy 


TTS 


2 poo 


(24) 


= I dr — 2 I <j{y)e ^ dy . 


'0 Jo 

By dominated convergence and due to (122|) the integral 2 a{y)e~^'' dy converges to some 

positive number in the limit £ —)■ 0. If we can show that the hrst integral converges to 0, 
the statement follows since then q[(p£] < 0 for some £ small enough and (Tess(—Ao-) = [0, cxd) 
by Theorem IIII.ll 

Regarding the hrst integral we hrst employ an iterated integration by parts to yield 


on—2 poo 

dr = --- / dr 


,2£-l„-2r« 

t \ I' ^ 

10 {n-iy.Jo 


(25) 


for any n > 2. Furthermore, employing a change of variables {y = r^) and setting e -.= n ^ 
we obtain 


POO POO 

I dr = n J e~‘^^y^~^ dy 


(26) 


= n2-^T{n) . 

Consequently, writing £ := n~^ and combining fl25|) and (|26|) . we can now calculate the hrst 
integral in (|2^ : we obtain 

-2 PCX) ^ on—2 /-OO 

lo 


I 2^—^ —2r^ i TT 2 

^26 Ig 2r ^ ^ . 

2 .In 2n^n-l)\Jo 


e dr 


TT 

8n 


(27) 


This shows that the hrst integral in 0241) converges to zero as n —)■ oo, hence proving the 
statement. □ 


In a next step, assuming that supp a{y) C [0, L] for some L > 0, we estimate the number 
of eigenstates of — below the essential spectrum. As in the proof of Theorem lIII.il we use 
a comparison argument between —A„ and now a direct sum of Laplacians —Aq^ © —Aq^ 
where = ffi U and where ffi is the triangle which is obtained by connecting the 
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three points (0,0),(L,0) and (0,L). On an operator level, we impose Nenmann boundary 
conditions along the line connecting (L, 0) and (0, L) as well as on 8 ^ 2 - Furthermore, Robin 
boundary conditions as in (jl]) with constant a := ||o'(?/)||oo are imposed along the other sides 
of the triangle Oi. 

On a form level, the quadratic form associated with —is given by 

:= / \V(p\‘^dx-a[ \ipbv{y)\^ dy (28) 

Jfli JdRlndfli 

and it is defined on Furthermore, the quadratic form associated with — is given 

by 

94 ^ := [ \Vip\^ dx , (29) 

J ^2 

being defined on H^(Q 2 )- As a consequence, the operator — Aq^ © — Aqj is (in the sense 
of an operator bracketing) smaller than Hence, denoting by N_(A) the number of 

negative (discrete) eigenvalues of a self-adjoint operator A, we obtain the following result. 

Lemma III.4. Let —A„ be given with supp <j(y) C [0, L] and corresponding — A^^ as con¬ 
structed above. Then 

N_{-A^) < N_(-An,) . (30) 


Proof. As mentioned already, in the sense of an operator bracketing we have —An, 0—A p., < 
—Ao-. A direct application of the minimax principle (see, e.g.. Corollary 12.3 in Ref. IlSh 
hence implies N_(—A„) < A_(—Ap^ © —Ap^). The statement then follows by taking into 
account that A_(—Ap^ ©—Ap^) = A_(—ApJ, noting that —Ap 2 is a positive operator. □ 


There is now an interesting way to estimate A_(—ApJ further by reducing the two- 
particle (or two-dimensional) problem to a one-particle and hence one-dimensional problem. 
For this, let — A^- denote the one-dimensional Laplacian — ^ being defined on 

V(—A.o-) = G H^(0, L) I ‘pfd) + dg)(d) = 0 and — ‘p'(L) + dip(L) = 0 } , (31) 


with (j > 0 (this operator is a “compact” version of the operator considered in Remark llll.2j) . 


It is well-known that (—A^-,!) 
negative energy (see, e.g.. Ref. 
corresponds to the solution of 


f-As)) is self-adjoint and has at most two eigenstates at 
161) . Furthermore, the ground state eigenvalue Sq := —k^ 


= a 


Ktanh 


k,L 

~ 2 ' 


10 


(32) 









with K > a > m. 

The next step is then to extend the triangle hli on which the operator — is defined in 
a suitable way: i.e., hli is extended by refiection across the line connecting (L, 0) and (0, L) 
to a square hli of side length L. By construction, Robin-boundary conditions with constant 
a are imposed along the sides of hli. However, the Neumann boundary conditions which 
were present along line connecting (L, 0) and (0, L) are now implicitly implemented through 
requiring symmetry of the functions across this line. To repeat and to be more precise, we 
define the operator 

- 0 1 + 1 0 -A^ (33) 

on the Hilbert space where hi = [0,L] x [0, L] and where the index s refers to the 

fact that only functions which are symmetric across line connecting (L, 0) and (0, L) are 
considered, i.e., f{x,y) = f{L — y,L — x). By symmetry it is readily verified that the normal 
derivative of any (differentiable) function / vanishes along the diagonal y = L — x. 

If a{—Aa) = {Sn I n G No} denotes the spectrum of —A^- and {SnlngNo the corre¬ 
sponding eigenvalues in increasing order and counted with multiplicity, then 


(t(— As- 0 1 -f 1 0 —As-) = {sn + Em \ n,m E No and n > m} 


(34) 


and {Sn + Em I n,m E No and n > m} are the eigenvalues of —As- 0 1-1-10 —As-, again 
counted with multiplicity: indeed, if ^-^e the eigenfunctions of —As-, each (tensor) 

product forms an eigenfunction of —As- 01-1-10— As- to the eigenvalue En+E^ if the 

underlying Hilbert space is L‘^{Q) (see, e.g.. Ref. 1 171). However, since we are only working on 
the symmetric subspace we only have to consider symmetrised eigenfunctions, i.e., 

eigenfunctions of the form ipn 0 ipm + ® ^n- As a consequence, for each n,m E Nq, the 

two pairs ipn 0 y^m and ipm ® y>n yield only one symmetric eigenfunction and, in order to 
avoid an overcounting of eigenvalues, one therefore has to require n > m. 

Since Eq = —tA we obtain, assuming that (— As-, P(— As-)) has exactly one eigenstate at 
negative energy, 

A_(-AqJ < #{n G No I ■ (35) 

Furthermore, employing the results obtained in Ref. [la, we see that (— As, T’(— As-)) has 
only one eigenstate at negative energy if and only if 

2 


a < 


(36) 
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Hence, summarising the conclusions from above we have established the subsequent state¬ 
ment. 


Theorem III.5. Let — A^. be given with supp(T(?/) C [0,L] and a = ||cr(|/)||oo < If 
denotes the one-dimensional Laplacian as constructed above with eigenvalues {£,i}neNo; one 
has 


iV-(-A,) < #{ii e No I £„< ko|} . 


(37) 


It is interesting to note that, using equation (4.4) of Ref. 


19 


, all positive eigenvalues of 


(—Aa-, T>(—Act)) can be written as Sn = k'^ > 0 where kn is a (positive) solution of 


tan(fcL) 


2ak 


(38) 


As a consequence, combining equations fl3^ and fl38|) and using Theorem IIII.5I we arrive at 
the following. 


Corollary III.6. For any fixed value L > 0 there exists (Tq > 0 such that —A^ has at most 
one eigenstate at negative energy for potentials a{y) with supp(j(?/) C [0,L] and ||a(j/)||oo < 
(To. If, in addition, <j{y) fulfils (|22|) then there exists exactly one eigenstate. 


C. On the ground state energy 

In this hnal subsection we provide upper bounds on the bottom of the spectrum. This 
result then provides us, in particular, with an estimate of the lowest eigenvalue below the 
essential spectrum in the case where such an eigenvalue exists. 

Theorem III.7. Assume cr{y) G L°°(R+) and set a := ||(t(?/)||oo- If Ea = inf (t(— Act) 
denotes the ground state energy then 

poo 

- 2 (7^ <E^< -2a^ + 8a^ / [a- cr( 2 /)]e“^"^ dy . (39) 

Jo 

Proof. The lower bound follows from a comparison argument: For this consider the operator 
—Aa- with a{y) := a as described in Remark IlII. 21 We readily see that this operator is indeed 
smaller (in the sense of an operator bracktering) than the given operator —Act with aly). 
Furthermore, since the lowest spectral value of — Aa is the eigenvalue —2d^, the lower bound 
follows directly. 
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To obtain the upper bound we use the (normalised) test-function ip{x,y) = 2ae 
i.e., the ground state of the example discussed in Remark IIII. 21 We calculate 


dy 


g[(p] = 2(j^ — 8(T^ / (y{y)e 

Jo 

poo 

= —2(j^ -I- 8(7^ / [a — cr(j/)]e“^°'^ dy . 

Jo 


This proves the claim. 


( 40 ) 


□ 


(^{y) ■= 


(41) 


We can illustrate the result of Theorem IIII. 71 by the following example: Consider the 
step-potential 

a ii 0 < y < L 
0 ii y > L 

where a > 0 is the depth of the potential and L > 0 the range of it. In this case we know 
that the essential spectrum starts at zero and that there exists a negative eigenvalue (see 
Theorem IIII. II and Theorem IIII. 3^ . Consequently, Theorem IIII. 71 implies 


2a^ <E„< -2a^ + 


(42) 


leaving us with an estimate of the lowest eigenvalue of the system. 


IV. ASYMPTOTICS OF THE GROUND STATE 


In this hnal section we will discuss certain properties of the ground state (po{x) (where 
X G M+) of the system. Note that ipo G 77^(M^) with ||</3 o||l2(r^) = 1 is called ground state 
of the system if 

g[(po] = E„ = inf cr(-A<^) > -oo . (43) 


In particular, we want to describe its asymptotic behavior as |x| —)■ cxo in the case where the 
boundary potential a{y) has compact support. 

Adopting standard methods available from the theory of Schrodinger operators and, in 
particular, the methods used in the proof of Theorem 11.8 in Ref. 
the following result. 


20 


, we directly arrive at 


Lemma IV. 1. Assume that ipo{x) is a ground state of Then (po{x) is unique (up to a 

constant phase) and can he chosen to he real-valued and strictly positive in Furthermore, 
(po{x) fulfils the eigenvalue equation — Ao-</3o = E^ip^. 
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As a final result we derive (pointwise) upper bounds on ipo{x) that describe its asymptotic 
behavior as |a;| oo assuming, for simplicity, that (T{y) has compact support. The key 
ingredient is an approach based on the concept of (positive) supersolutions as described in 


Ref. 


21 


More precisely, if < 0, the supersolution in our case will be the function 

^-y/\E7\\x\ 

u{x) = 




(44) 


considered on Qpi := {x G | |a;| > R} for some R > 0 large enough. It is called a 
supersolution since it fulhls (—A + > 0 on fl/j. In the proof we will also use the 

fact that the eigenfunction ipo can be extended (by reflection across the coordinate axes) to 
all of Qr and assumed to be continuous. Indeed, for R large enough, the boundary conditions 
on are Neumann boundary conditions only (here enters the assumption of compact 

support). Hence, standard regularity theory (see, e.g.. Ref. l22|) implies that is 

in n Hh)- After having extended (fo to by reflection and then, using a cut-off 

function, to all of standard Sobolev embedding theorems imply continuity (see, e.g.. 


Ref. 


23|). Note that extending ipo to all of while remaining in fl Qr) is possible 


due to the Neumann boundary conditions (vanishing normal derivative). 

Theorem IV. 2. Assume that supp(j(j/) C [0,L] for some L > 0. If (po is the ground state 
of — Ao- with ground state energy E„ < 0 then, for each R> L, there exists a constant c > 0 
such that 


|<y9o(x)| < C 


,-R\E,,\\x\ 




Vx e 


|x| > R -I-1 


(45) 


Proof. We only sketch the proof. For more details see Theorem 2.7 and Theorem 3.2 in 
Ref. I 21 I As described above, we will assume (po{x) to be extended to and continuous on all 
of Qr. Now, picking R > L and dehning Rq := R + 1 we choose c > 0 such that 


cm(x) — (po{x) >0, Vx : |x| = Rq • (46) 

Note that (IT6]) is meaningful due to continuity of both functions u{x) and (po{x) on {x | |x| = 
Ro}. The key idea now is to prove that the function 

uo{x) := [y:>o{x) -cn(x)]_^ (47) 

vanishes on all of VLr^. Note that [•]_!_ denotes the positive part of a function. 
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From continuity we conclude that there exists 5 > 0 such that uo(x) = 0 for all x 
s.t. Rq < |a;| < Ro + S. Furthermore, we note that (po(x) — cu(x) is a subsolution, i.e., 
(-A + |E,|)(v„ (x) — cu{x)) < 0 on Qrq. Lemma 2.9 of Ref. l2l| then shows that Uo{x) is 
also a non-negative subsolution (in the weak sense). Hence, 


VmoV(C^mo) dx + \Eo 


dx < 0 , 


(48) 


‘Ro 


‘Ro 


for any real function ( G C^(f2/{Q). Using the identity (see equation (2.25) in Ref. l2ll) 


VtioVK^uo) = |V(Ctio)t - ti?|VC|- 


(49) 


we obtain 


(|V(Cmo)P + l^alC^Mo) da: < f MqIVCP dx 

^Rq - 


Now, using fl5U]) and the identity (see equation (2.6) in Ref. l2l[l 


VnV 


^)=|V(C«o)t-«^ 

u 


Cuo 


u 


we obtain 


u 


'■Ro 


-(^) 


dx < / dx , 


in 


Rq 


(50) 


(51) 


(52) 


in analogy to equation (2.28) in Ref. l21. 

We then choose a suitable sequence (Xn)nGN with x„(x) := x(^), X ^ C'“(M^), 0 < y < 1 
and x(x) = 1 for |x| < 1. Employing the Lemma of Fatou as well as the estimate (|5^ we 
obtain 

2 


U 


'^Rq 


) dx<liminf / dx 

\uJ n^oo 


(53) 


= 0 . 


As a consequence, Mo(a:) = \u{x) with some constant A. Since u > 0 in while uq 
vanishes in some strip around Rq it follows that A = 0. Thus Mo(a:) = 0 and the theorem is 
proved. □ 

Remark IV. 3. Theorem IIV.2l establishes the asymptotics of the ground state wave function 
(po{x) as |x| oo, assuming the boundary potential a{y) has compact support. In order 
to understand how the compactness of the support affects the asymptotics, it is instructive 
to look at the ground state of the example discussed in Remark 1111.21 i.e., the case of 
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aiy) = a > 0. In this case the ground state is given by (p{x,y) = Setting 

X = y hr simplicity, one has \ip{x,x)\ = taking into account that 

|a;| = \/2x and that E„ = —2a‘^ is the corresponding ground state energy. We hence conclude 
that the finite range of the boundary potential leads, in the direction where x = y, to an 
additional decay factor of 1/ \/\x\. 
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